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ABSTRACT 

Spectral absorption features in active galactic nuclei (AGNs) have traditionally been attributed to outflowing 
photoionized gas located at a distance of order a parsec from the central continuum source. However, recent 
observations of QSO FIRST J104459. 6+365605 by de Kool and coworkers, when intepreted in the context of a 
single-phase gas model, imply that the absorption occurs much farther (« 700 pc) from the center We reinterpret 
these observations in terms of a shielded, multiphase gas, which we represent as a continuous low-density wind 
with embedded high-density clouds. Our model satisfies all the observational constraints with an absorbing gas 
that extends only out to 4 pc from the central source. The different density components in this model coexist 
in the same region of space and have similar velocities, which makes it possible to account for the detection in 
this source of absorption features that correspond to different ionization parameters but have a similar velocity 
structure. This model also implies that only a small fraction of the gas along the line of sight to the center is 
outflowing at the observed speeds and that the clouds are dusty whereas the uniform gas component is dust free. 
We suggest that a similar picture may apply to other sources and discuss additional possible clues to the existence 
of multiphase outflows in AGNs. 

Subject headings: galaxies: active — galaxies: Seyfert — quasars: absorption lines — quasars: individual 
(FIRST J104459.6+365605) — MHD 



1. INTRODUCTION 

I Modeling active galactic nuclei (AGNs) is a challenging en- 
■ deavor in part because of the uncertain geometry of the gas 
' and dust within their cores. Our best understanding of the 
gas distribution comes from reverberation mappings (Bland- 
ford & McKee 1982; Netzer & Peterson 1997), which indicate 
I that the broad emission-line region (BELR) lies at a distance 
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' ^BELR ~ O.OIL44 pc from the central continuum, where L44 is 
, the luminosity in the spectral band ^0.1-1 fim in units of lO'*'* 
' ergs s~'. Based on observed absorption of BELR features by 
broad absorption-line region (BALR) gas and on optical polar- 
ization measurements, it has been inferred that the BALR lies 
outside of the BELR; the best estimates have placed the BALR 
' somewhere between a few tenths and several parsecs from the 
central source (e.g., Turnshek 1988). 

Velocities for these outflows are well determined from 
Doppler shifts: the BELR and BALR components have speeds 
< 5000 to 8000 km s"' and < 30000 km s"' , respectively. Other 
outflow components have also been detected. Observations 
of the "warm absorber" (a partially ionized X-ray absorbing 
gas) in Seyfert galaxies and radio-quiet QSOs, and of a UV- 
absorbing component apparently associated with the warm ab- 
sorber, indicate gas outflowing at < 10^ km s~'. Furthermore, 
spectral lines classified as being from the Narrow Line Region 
(NLR) correspond to velocities of < several hundred km s"' . 

Different theories have been proposed to explain some of 
these outflows. For example, Emmering, Blandford, & Shlos- 
man (1992) considered a hydromagnetically-driven disk out- 
flow of discrete, long-lived clouds as a model of the BELR, 
whereas Murray et al. (1995) and Proga, Stone, & Kallman 



(2000) explored radiatively-driven, continuous disk winds as 
the origin of the BELR and the BALR. These models all in- 
volved a single-phase gas medium. 

Recently, de Kool et al. (2001, hereafter dKOl) analyzed 
a Keck HIRES spectrum of the FIRST Bright Quasar Source 
J104459.6+365605 (hereafter, FBQS 1044). They inferred that 
excited Fe II levels are not populated according to local ther- 
modynamic equilibrium. This observation implies an electron 
number density n^, « 4 x lO^* cm"-', which, together with the 
observed broad absorption in Mg II [qualifying this object as a 
broad absorption-line quasar (BALQSO)], the Mg I absorption 
features, and single gas-phase photoionization models, place 
this BALR at ^ 700 parsecs from the central source. 

Not only does this surprising result contradict the picture of a 
close-in BALR, but, in addition, dKOl's thorough analysis sug- 
gests that the gas only partially covers the continuum source. It 
is not easy to understand how the absorbing gas could be ap- 
proximately 1 kpc from the center and yet only partially cover 
the AGN core. Additionally, the Keck spectra indicate distinct 
groupings of gas in velocity space; it is also difficult to ex- 
plain how these apparent clumps could move to such a large 
distance and still remain bunched up, retaining distinct identi- 
ties. Even more striking is dKOl's observation that different 
ionization states have a similar kinematic structure (i.e., resid- 
ual intensity as a function of velocity): in particular, both Fe II 
and Mg I exhibit absorption troughs at -200, -1250, -3550, and 
-3800 km s-i. 

The 700 pc distance estimate results from combining the ion- 
ization parameter^ implied by the observed low-ionization Mg 
and Fe absorption lines with the inferred values of rtg and L 
(= lO'^'ergs s-i). One can reduce this distance estimate by atten- 
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uating the incident continuum by an intermediate gas "shield." 
In this model, however, ions at different ionization states (such 
as the Fe II and Mg I components seen in FBQS 1044) are ex- 
pected to occupy different regions in space, with lower val- 
ues of U corresponding to larger distances from the center 
(e.g., Voit, Weymann, & Korista 1993). This stratification will 
likely lead to disparate velocities for the different ions: for ex- 
ample, models of self- similar MHD winds and of radiatively- 
driven constant-t/ outflows predict that the terminal velocity, 
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Voo, varies as '"inject) where inject is the radius at which the gas is 
injected into the outflow (e.g., Blandford & Payne 1982; Arav, 
Li, & Begelman 1994). Models that rely solely on shielding and 
distance to separate the different ionization components there- 
fore cannot explain dKO 1 ' s observations, in which different ion- 
ization states are found to have similar velocities. 

We propose to explain the observations of FBQS 1044 by 
generahzing the single-phase shielded-gas model, attributing 
the different ionization states to different density components 
in a multiphase outflow. These components coexist at the same 
distance from the center and thus have different ionization pa- 
rameters but essentially the same velocities. In particular, if the 
high- and low-ionization lines arise, respectively, in a continu- 
ous wind and in dense clouds that are embedded in the outflow, 
then all the absorption components produced in a given region 
of the wind will exhibit similar kinematic signatures. 

In our model, the continuous gas component extends from 
near the black hole's event horizon out to the distance where 
the observed Fe 11 absorption and electron density can be re- 
produced. The inner part of this component (interior to the 
BALR) is identified as the "shield." This region could be asso- 
ciated with an MHD-driven (Konigl & Kartje 1994), a Thom- 
son scattering-driven (e.g., Blandford 2001), or a "failed" hne- 
driven (e.g., Murray et al. 1995; Proga, Stone, & Kallman 2000) 
disk wind, or with a disk corona (e.g., Emmering et al. 1992). 
The outer region of the continuous gas component is outflow- 
ing and accounts for the Fe II and Mg 11 absorption, but is still 
too highly ionized to contain Mg I absorbing gas. Within this 
outflow are embedded higher-density clouds that account for 
the lower-ionization Mg I absorption: their high density yields 
a lower U in the clouds, allowing Mg I to exist. Such a two- 
component outflow may arise naturally in the context of a cen- 
trifugaUy driven disk wind, which could uplift clouds from the 
disk surface by its ram pressure and confine them by its inter- 
nal magnetic pressure (e.g., Emmering et al. 1992; Kartje et 
al. 1999; Everett, Konigl, & Kartje 2001). Alternatively, the 
clouds may represent transient density enhancements that are 
produced by turbulence (e.g., Bottorff & Ferland 2001) or by 
shocks in a radiatively driven wind (e.g., Arav et al. 1994). For 
illustration, we adopt here the "clouds upUfted and confined by 
an MHD wind" picture. 

2. MULTIPHASE OUTFLOW MODEL 

2.1. Setup 

Our model consists of two segments that represent two dis- 
tinct gas phases. The first segment, which corresponds to the 
continuous phase, extends from rin = 10 GM/c^ = 1 A x 10'^ cm 
(following dKOl, we take the central black-hole mass to be 
M = 5 X Kf'Mo) out to just past the hydrogen recombina- 
tion front, which corresponds to the Fe III Fe II transi- 
tion. The second segment, which models the confined dense 
clouds, is a constant-density zone located just beyond the con- 
tinuous gas component (since the cloud temperature is nearly 



uniform, this provides a good representation of magnetically 
confined constant-pressure clouds). In reality, the continuous 
outflow in the line-absorption region can be expected to engulf 
the clouds, but this simplified model should capture the basic 
physical effects of a two-phase medium. (We verified that our 
results are independent of the cloud radial distribution within 
the Fe II/Mg I absorption region.) 

We obtain the ionization structure of the different gas phases 
using the photoionization code Cloudy (Ferland 2000), adopt- 
ing the same source parameters as in dKOl and assuming (as 
they do in most of their models) a Mathews & Ferland (1987) 
spectrum. We take the scaling of the hydrogen number den- 
sity with spherical radius to be n(r) oc r"', as in the AGN disk- 
wind models considered by Konigl & Kartje (1994). In these 
self-similar MHD outflows, the magnetic field Bwind also scales 
as r~\ which results in magnetically confined clouds having a 
constant ionization parameter. 

We first calculate the photoionization of the continuous com- 
ponent, stopping the computation when tie drops to 4 x 10^ 
cm"^ (the value inferred from the observations). We then cal- 
culate how the radiation emerging through the continuous seg- 
ment affects clouds located at that distance. The two parameters 
we adjust in our model are n^j, the wind hydrogen number den- 
sity at Tin, and ric, the hydrogen number density of the clouds. 
We explore a range of models to find the values of and ric 
that best reproduce the observations. 

2.2. Results 

Our best-fit value for the wind density is n„^i « 10**'^^ cm"^. 
This yields the observed tie in the region {al r k, A pc) where 
the Fe II column density in the wind attains the observed value 
of 3 x 10^^ cm"-^, which is significant since the measurement 
comes from the Fe II absorption lines. However, we do have to 
cut off the outflow very close to the end of the hydrogen recom- 
bination front so as not to exceed the observed Fe II column 
(see Fig. 1). Although the abrupt end of the absorbing colunm 
could be an artifact caused by our simplified treatment, the oc- 
currence of a strong reduction in the absorption at this location 
may have an actual physical basis. It may be a consequence 
of the decrease in the ionization fraction in this region, which 
reduces the efficiency with which the disk can drive an MHD 
outflow, and it may also reflect a transition from a gaseous to a 
"clumpy" disk (see Shlosman & Begelman 1987) or from dif- 
ferential to sohd-body disk rotation (see Bottorff, Korista, & 
Shlosman 2000) on that scale. 

We also note from Figure 1 that the Mg II front occurs ~ 2.5 
times closer to the central source than the Fe II front. (This is 
due to the formation of a He III He II front in the wind that 
absorbs photons with E > 54.4 eV, allowing Mg III, with an 
ionization energy of 80.1 eV, to recombine into Mg II.) Since 
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Voo OC '"inject i'^ °™ choscu disk-wiud model (see §1), the pre- 
dicted smaller initial radius of the Mg II front is consistent with 
the higher Mg 11 velocities observed in FBQS 1044. 

Besides determining the measured values of the Fe II col- 
umn and of tie, the continuous wind component in our model 
also exerts a strong influence on the cloud absorption proper- 
ties through its effect on the transmitted continuum spectrum. 
We find that, in order for the clouds to account for the observed 
Mg I absorption without also dominating the Fe II absorption, 
their iron gas-phase abundance must be reduced, which we at- 
tribute to depletion by dust. We assume an ISM dust compo- 
sition [the default Cloudy graphite and siUcate grains with a 
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Mathis, Rumpl, & Nordsieck (1977; hereafter MRN) power- 
law size distribution]. The dust abundance is constrained by 
the measured extinction: A2500 ^ 1 - O- 

Our best- fit value for the typical cloud density is ric « 10^-^ 
cm"^. This means that the clouds are ^ 10"* times denser (and 
have a correspondingly lower ionization parameter) than the 
surrounding continuous wind. Figure 2 shows the absorption 
column densities associated with a string of such clouds as a 
function of pathlength through the clouds. 

Our model can account for the observations with cloud den- 
sities that range between 10^^^ and 10'^ cm"^. This leeway re- 
flects the uncertainty in A25oo- Whereas dKOl employ an ex- 
tinction of 1.0 magnitude, A2500 could range from ^1.0 down 
to 0.2 magnitudes based on a comparison between the spectrum 
of FBQS 1044 and that of an "average" QSO (Weymann et al. 
1991). If a population of lower-density clouds were present 
with a significant fraction of the observed Mg I colunnn, then 
the dusty clouds would have A25()o > 1 . On the other hand, if 
He > lO** cm"^, A2500 would drop below 0.25. With ric = 10^-^ 
cm"^ we predict A2500 ~ 0.5 magnitudes, which fits weU within 
the inferred range. 

With the extinction already accounted for, the data require 
the continuous gas component to be effectively dust free. We 
tested this hypothesis by including Cloudy's Orion-type dust (a 
set of graphite and silicate grains with MRN's minimum size in- 
creased from 0.0025 to 0.03 /xm, appropriate to a UV-irradiated 
medium) in the wind, taking into account dust sublimation and 
sputtering. To include the effects of the latter process, we re- 
moved all dust grains whose sputtering time scales (calculated 
following Tielens et al. 1994) do not exceed 1000 yr (the ap- 
proximate time grains spend in the wind if they travel at the 
observed outflow velocity over a distance of 1 pc). Including 
both of these effects, the predicted dust extinction is over one 
order of magnitude greater than the maximum value allowed 
by the observations. We also considered the effect of chang- 
ing the scaling of the wind density with radius from n oc r"' to 
n (X r"^/^ (as in Blandford & Payne 1982) and to n oc r~^. These 
models again overestimate the extinction by about an order of 
magnitude after dust subUmation and sputtering are taken into 
account. We therefore conclude that, if the gas comprising the 
shield is associated with a disk outflow, then it must already be 
dust free when it leaves the disk. Such a situation could arise 
if the wind originates in a hot disk corona where the matter re- 
sides long enough for any dust grains to be destroyed. 

However, to provide the requisite shielding, a dust-free wind 
requires a higher total gas column than a dusty outflow, so the 
inferred lack of dust in the wind implies a prohibitively large 
mass outflow rate — much larger than could be launched by 
the magnetic field that confines the clouds in our model. Thus, 
only part of the shield can be outflowing at velocities that are 
comparable to (or exceed) the observed speeds. Perhaps the in- 
ner outflow has a lower velocity than we predict, or else some 
of the shield may not even be outflowing, as in a disk corona or 
a "failed" line-driven wind (see §1). 

All of these considerations yield our best model, which is 
compared with the observational results in Table 1 . This model 
satisfies aU the observational constraints and impUes that the 
absorbing gas lies over two orders of magnitude closer to the 
central source than the earlier estimate. The radius where the 
observed gas leaves the disk surface is, in general, smaller yet. 

We can estimate the mass outflow rate associated with the 
absorbing gas through the relation Mwind ~ '^i^rfNyimpV, where 



/ is the fraction of 47r steradians into which the wind flows, 
A^H is the total hydrogen column density of only the inferred 
Mg II and Fe II absorbing region (w 3.9 x 10^^ cm~^), v is 
the observed outflow speed (~ 10^ cm s~'), r (= 1.2 x 10^^ 
cm) is the inferred distance, and where we assume a verti- 
cally and azimuthally continuous wind. We take / ~ 0.1 for 
BALQSO sources (Weymann 1997). For the above values, we 
find Mwind « 8Mq yr"i. 

By equating the thermal pressure in the clouds (« 2.7 x 10""* 
dynes cm"^) to the confining wind magnetic pressure, B\^^^^/%it, 
we deduce 5wmd ~ 8.2 x 10"^ G. It is encouraging that, when 
this value of Bwind is used in the n{r) oc self-similar MHD 
wind model, it implies a local mass outflow rate that is com- 
parable to the above estimate of Mwind, yielding ^ 4 Mq yr"^ . 
If we instead choose = 10^ cm"^ (so as to satisfy the lower 
limit on the dust extinction, A2500 = 0.25) and require pressure 
balance, we find Mwind ~ 9 Mq yr~^ 

3. CONCLUSION 

The most robust result of our study is that a shielded, multi- 
phase absorption region reproduces the observations of FBQS 
1044 on a conventional BALR scale (« 4 pc). In addition, when 
one attributes the Fe 11 and Mg 11 absorption to a low-density 
outflow component and the Mg I absorption to a cospatial high- 
density outflow component, it is possible to explain the simi- 
lar kinematic structure of the respective spectral features. We 
also find that only a small fraction of the gas along the line of 
sight can be outflowing at the observed speeds and that only 
the high-density component of the outflow is dusty. We derived 
these results using a "clouds embedded in a continuous MHD 
disk wind" model, but our conclusions also apply to other plau- 
sible scenarios that include a continuum-shielding gas column 
and an absorption region that contains distinct low- and high- 
density components. Our basic conclusions appear to be quite 
general, although the precise composition of the absorbing gas 
and its detailed spatial and kinematic properties are not fully 
constrained by the observations and remain model dependent. 

In addition to explaining the FBQS 1044 observations, this 
picture may be relevant to the interpretation of absorption fea- 
tures in similar objects where single-phase models imply a large 
distance. For instance, in the case of the radio-loud galaxy 3C 
191, absorber distances of ~ 28 kpc were inferred by Hamann 
et al. (2001) using similar arguments to those employed by 
dKOl. A multiple-phase model could place these absorbers 
much closer to the central source. 

This interpretation may also be apphcable to other AGN ob- 
servations. As outlined in § 1, several distinct outflow com- 
ponents have been inferred in various types of AGNs. There 
is now growing evidence that these components may not be 
single-phase. For the warm-absorber component, which has 
been inferred to give rise to both X-ray and UV absorption 
(e.g., Crenshaw 1997; Mathur et al. 1998; Monier et al. 2001), 
there are indications in at least some sources that the X-ray and 
UV absorbing components are not identical (e.g., the Seyfert 1 
galaxy NGC 3783 — Kaspi et al. 2000; Kraemer, Crenshaw, 
& Gabel 2001). Furthermore, it appears that the UV and X- 
ray absorbing components are themselves divided into multiple 
zones. In the UV regime, this has been established in Seyfert 1 
galaxies like NGC 3783 (Kraemer, et al. 2001) and NGC 3227 
(Crenshaw et al. 2001). In the X-ray regime, this is exempli- 
fied by the Seyfert 1 galaxy MCG-6-30-15, which was mod- 
eled by Morales, Fabian, & Reynolds (2000) as an extended 
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multi-zone medium. These authors suggested that, in reahty, 
these zones may correspond to a continuum of clouds at dif- 
ferent radii and different densities, as in the BELR model of 
Baldwin et al. (1995). A hierarchical, turbulent-gas realization 
of this concept was recently presented by Bottorff & Ferland 
(2001). Our "clouds in a continuous wind" scenario is designed 
to explicitly address the outflowing nature of the absorbing gas 
but is otherwise similar to the above picture in its description 
of spatially coexisting multiple phases. For simplicity, we have 
treated the clouds as long-lived, pressure-confined entities, but 
it is entirely conceivable that the clouds are transient features 
that arise in a turbulent outflow. The same basic picture of a 
multiphase medium may thus apply to the gas in the BELR, the 
BALR (e.g., Arav et al. 1999), the warm absorber, and even the 
NLR (e.g., Komossa 2001). 

In conclusion, we have argued that the observations of FBQS 
1044 can be interpreted in the context of the standard BALR 
picture in terms of a gas outflow that consists of (at least) two 



phases, which we modeled as a continuous wind with embed- 
ded dense clouds, shielded from the central continuum. If the 
shield is identified with the continuous component, then only a 
fraction of it can be outflowing at a speed that approaches (or 
exceeds) the value in the Fe/Mg absorption region. We also 
deduced that the clouds are dusty but that the shield is effec- 
tively dust free. As far as we are aware, this is the first instance 
of a BALR outflow in which the data provide direct evidence 
for the existence of a multiphase medium. Together with other 
pieces of evidence, this result lends support to the view that all 
the major outflow components in AGNs may contain multiple 
phases. 

We gratefully acknowledge the initial work of and continued 
discussions with Martijn de Kool as well as helpful comments 
from Mike Brotherton, Abraham Loeb, and the referee. I.E. 
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Table 1 

Comparison of the Two-Phase Outflow Model with Observational Results for FIRST QSO 1044-1-3656 



Quantity 


Observational Result 


Model Result 


N(Fell) 


~ 3 X lO'^cm"^ 


3 X lO'^cm"^ 


N(Fel) 


< lO'^cm-2 


2.5 X 10"cm-2 


N(Mg II)/N(Mg I) 


^30 


500 


N(Mgl) 


~ 2 X lO'^cm-2 


2 X lO'^cm-^ 


He 


~4000 


4000 


A25OO 


<1 


~0.5 
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Fig. 1 . — Column densities in tlie dust-free continuous wind model parameterized by n,,., ft; 10^^' cm"^. The continuous wind accounts for the Fe n as well as 
the higher-velocity Mg 11 absorption without overproducing the low-ionization Mg I absorption (which is attributed to the dense clouds). 
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Fig. 2. — Column densities in dusty clouds of hydrogen density Uc = 10*'^ cm ' that are shielded by the continuous dust-free wind of Fig. 1. The clouds account 
for the Mg I absorption without overproducing the higher-ionization Fe 11 absorption (which is attributed to the continuous wind). 



